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Abstract. Light traveling through a liquid crystal with disclinations perceives a geometrical background 
which causes lensing effects similar to the ones predicted for cosmic objects like global monopoles and 
cosmic strings. In this article we explore the effective geometry as perceived by light in such media. The 
comparison between both systems suggests that experiments can be done in the laboratory to simulate 
optical properties, like gravitational lensing, of cosmic objects. 

PACS. 61.30.Jf Defects in liquid crystals - 61.30.Dk Continuum models and theories of liquid crystal 
structure - 78.20.Bh Theory, models, and numerical simulation - 02.40.-k Geometry, differential geometry, 
and topology 



1 Introduction 

Topological defects may be formed when a symmetry is 
broken in a phase transition. Examples run from Cosmol- 
ogy to Condensed Matter Physics. In Cosmology, defects 
like global monopoles and cosmic strings might have ap- 
peared at the beginning of the universe when the sym- 
metries of the unified forces of nature were broken in a 
sequence of phase transitions while the universe cooled 
down. In Condensed Matter Physics, the prototype topo- 
logical defect is the Abrikosov-Nielsen-Olesen vortex 2 in 
type II superconductors with broken U(l) symmetry. On 
the other hand, uniaxial nematic liquid crystals, in their 
isotropic phase, have global SO (3) rotational symmetry, 
which is broken under cooling, yielding the nematic phase 
and producing topological defects like the hedgehog point 
defect or disclinations as relics of the former, less ordered, 
phase. Remarkably, they are formed according to the same 
physical principle, the Kibble mechanism [3], as the cos- 
mic defects, as verified experimentally by Bowick et al. 

H 

The similarity between defects in nematics and cosmic 
defects goes beyond their processes of formation. Cosmic 
strings, for instance, although do not introduce curvature 
(gravitational field) globally, they do change the topol- 
ogy and, therefore, introduce interesting phenomena like 
gravitational lensing 0. In this article, we use the ge- 
ometrical model developed by Joets and Ribotta [H] for 
anisotropic media to study the propagation of light in the 
nematic phase with topological defects, finding lensing ef- 
fects. This result, in fact, is not new. It has been previously 
observed as early as 1919 by Grandjean |Jj in a different 



context. Our aim in this work is to reveal the geomet- 
ric structure underlying this system and its similarity to 
cosmological defects, which permits that experiments may 
be performed with liquid crystals to study otherwise un- 
reachable objects like the global monopole or the cosmic 
string, in the spirit of reference 0|. 

Another system that has been compared with the cos- 
mic string is the irrotational hydrodynamical vortex [H]. 
Propagation of phonons in a medium with such vortices 
can be described by an effective Riemannian geometry 
which asymptotically agrees with that of the massless spin- 
ning cosmic string. The phonon rays are geodesies in this 
effective geometry and a converging lens behavior has been 
found [S]- 



2 Geometry 

In an inhomogeneous, isotropic, medium Fermat's princi- 
ple, stating that the optical length of a ray between two 
given points be the shortest one joining these points, is 
nicely interpreted as a geodesic in a Riemannian space 
[Sj. A generalization of Riemannian geometry, Riemann- 
Cartan geometry, has been a valuable tool to the study 
of topological defects both in Gravitation and in Con- 
densed Matter Physics ^U]- Yet another generalization 
of Riemannian geometry, Finsler geometry 1 1 1 j , has been 
used to model the liquid crystal as an optical medium 
[Hj. This geometry takes into account anisotropy by con- 
sidering that the line clement depends not only on the 
location but also on the local orientation of the line to be 
considered. This seems a very reasonable approach to light 
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propagation through anisotropic media since the local ori- 
entation of the molecules is incorporated by the geometry 
perceived by the light rays. Accordingly, Joets and Ribotta 
interpreted the light rays in an anisotropic medium as 
geodesies in a Finsler space. Although this "anisotropic" 
Riemannian geometry appear as a natural, powerful, tool 
to model such problem, ordinary Riemannian geometry is 
enough |12| for highly symmetric molecular configurations 
like the ones presented by topological defects in nemat- 
ics. For the cases studied here Riemannian geometry gives 
exactly the same answers as the more general Finslerian 
approach. Therefore, for simplicity, we present our calcu- 
lations in the Riemannian approach in the body of this 
article, leaving the Finslerian approach for the Appendix. 

It is important here to make a distinction between ordi- 
nary and extraordinary light rays propagating in uniaxial 
nematics: in the former the electric field is perpendicular 
to the director n and in the latter the electric field has a 
component along n. Differently from the ordinary ray, in 
an extraordinary ray, the direction of the Poynting vector 
S in general differs from the direction of the wavevector 
k. This implies that the energy velocity differs from the 
phase velocity, leading to two different refractive indices 
jjj: the ray index N r , associated to the energy velocity, 
and the phase index N p , associated to the phase velocity, 
given by 

N 2 = n 2 Q cos 2 (3 + n\ sin 2 j3, (1) 



where (3 = (n,S) and 



Nl = 



n 2 cos 2 7 + ri 2 , sin 2 7 ' 



(2) 



where 7 = (n, fe). In both equations n and n e are, re- 
spectively, the ordinary and extraordinary indices 

Among all possible paths between the generic points 
A and B, Fermat's principle for the extraordinary rays 
grants us that the path actually followed by the energy 
is the one that minimizes the integral (see, for example, 
section 3.3.3.2 of 



The meaning of this equation is the following: the line ele- 
ment of the optical path in Euclidean space with refractive 
properties is identified with the line element of an effective 
geometry characterized by gij . 

The equation for the rays is the geodesic equation in 
Riemannian space |14| . 



~~dP 



En 



dx 3 dx 
'■ dt dt 



= 0, 



(G) 



where t is a parameter along the geodesic and r l - k are the 
Riemannian connections, or Christoffcl symbols, given by 



dgkm , dg. 
dx 3 



+ 



dgjk 

dx k dx m 



(J) 



Among the many geometrical properties that derive 
from the metric, the Ricmann tensor and the curvature 
scalar are the most important because they give a feel- 
ing for how curvature is distributed in space. For the spe- 
cific expressions for the Riemann tensor and the curvature 
scalar in terms of the metric and Christoffel symbols, we 
refer the reader to any good textbook on Riemannian ge- 
ometry such as |14j . 

The effective metric gij felt by the extraordinary light 
propagating in the medium with ray index N r is obtained 
from Eq. JSJ. The angle (3 is determined from the specific 
configurations of the director field n as follows. If the curve 
R(£), where t is a parameter along the curve, represents 
the light trajectory, then 



T{t) 



dR 



(8) 



is the tangent vector to the curve at each position parametrized 
by t (see Fig. 1). It is a well known result of the differential 
geometry of curves (see for example section II B of |15|1 
that = 1 if t is the arc length, which is the case. 

Since ||n|| = l, we have that 



cos/3 = T ■ n. 



(9) 



N r dt. 



(3) 



Here, dt is the element of arc length along the path. Ordi- 
nary rays are not interesting because their refractive index 
is simply n , a constant. 

In Riemannian geometry the line element ds depends 
on the position coordinates x l of the point of the manifold 
under consideration 



ds 2 = gjjdx l dx 3 . 



(4) 



The geodesic joinning points A and B in such manifold 
is obtained by minimizing J ds, just like Fermat's princi- 
ple. This leads to a nice interpretation of the light paths 
as geodesies in an effective geometry Thus, we may 
identify 

N 2 dt 2 =^g ij dx i dx 3 ' . (5) 

i,3 




Fig. 1. The curve R(£) in three-dimensional space, its tangent 
vector T(£) and the director n. 
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We study the hedgehog point defect and wedge discli- 
nations of generic strength k and constant c. The hedge- 
hog director configuration is characterized by the unitary 
radial vector 

n = f (10) 

in spherical coordinates, whereas the disclination configu- 
rations in the plane z = const, are given by 



n = (cos ip, sin<y9,0), 
in Cartesian coordinates, where |13| 
ip = ktp + c, 



(11) 
(12) 



where <fi is the angular polar coordinate. Selected director 
configurations can be seen in Figs. 2-4. 

It is assumed that the optical axis lies along the aver- 
age molecular orientation. Also, assuming the line defects 
straight, their translational symmetry along their axis al- 
lows us to ignore the dimension along the defect and study 
a two-dimensional cross section of the medium, the plane 
z — const. A word of caution here: the k = ±1 line de- 
fects are known to be unstable with respect to a escape 
into the third dimension 12|- Therefore, they might be 
very difficult to be observed. In fact, recent results 
indicate that even in thin films the escaped texture is the 
stable configuration. One might think that director fields 
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Fig. 

k = 



2. Director field for the hedgehog point defect and the 
-1, c = ? disclination, respectively. 
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similar to the k — ±1 defects could be obtained with the 
help of electric or magnetic fields or even by surface an- 
choring schemes but, to the best of our knowledge, this 
has not been observed so far. So, the results concerning 
the k = ±1 defects are, in this respect, just a theoretical 
possibility. For experimental procedures the hedgehog and 
k = ±b line defects will be more interesting due to their 
known stability 



2.1 Hedgehog 

Since the hedgehog is spherically symmetric we use spher- 
ical coordinates and write R — rf, in the spherical coor- 
dinate basis {f, 9, <j>}. Thus, 



dR 



dr 
dl 1 



df 
"dl' 



(13) 



But, using the Euclidean line element in spherical coordi- 
nates, 

dl 2 = dr 2 + r 2 (d9 2 + sin 2 6d<p 2 ) (14) 

and Eq. JTSJ and, since ^ is orthogonal ^2] to f, Eq. © 

becomes 

cos/3 = r, (15) 

leading to 

sin /3 = r\Je 2 + sm 2 6»</> 2 , (16) 

where f=jg, 6 = f t and^=f. 

Eqs. JT5J, Unjl, with the help of Eq. JTJ, make Eqs. |@J 
and into 



ds z = N z dl z = n z Q dr z + n z e r z {d6 z + shi 8d<j) z ). (17) 
Now, if we rescale the radial coordinate by making 



p = n r the metric (|17f) becomes 
ds 2 = dp 2 + a 2 p 2 (d0 2 



(18) 



with a = n e /n a . Notice that, in the isotropic limit, n Q = 
n e , the Euclidean line element in spherical coordinates, 
Eq. (|14l) is recovered. The metric (I18|l is identical to the 
space section of the global monopole metric ^7] . 

The nonzero Riemann tensor components correspond- 
ing to the metric ljl7JI is given by 



R, 



n 2 e (n 2 e 



"o) 2 -2, 



-r sin 



(19) 
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3. Director field for the strength k - 
^ disclinations, respectively. 
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Fig. 4. Director field for the k 
disclinations, respectively. 



and k 



and the permutations 
RothOch- The curvature scalar is 



-R, 



-R 



R = 



2{n 2 e 



(20) 



Again, in the n = n e limit, the Euclidean space (null 
curvature everywhere) is recovered. 

In summary, the effective geometry perceived by light 
in the presence of a hedgehog point defect, as described by 
its metric, Riemann tensor and curvature scalar, coincides 
with the space part of the geometry associated to a global 
monopole [T7| . a topological defect of spacetime. 
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2.2 Disclinations 



3 Light Paths 



Since we consider light propagating in the transversal z = 
const, surface to the line defects, we write R = r cos <frx + 
r sin combining the polar coordinates r and with the 
Cartesian basis {x, y}. Following the same procedure as in 
the previous subsection, we get 

T = {f cos — r<fi sin 0)i + if sin + r0 cos 0)y, (21) 

where, again, r = =jj and </ > = Now, with the help of 
Eqs. fl"T|) and Eq. © becomes 

cos j3 = r cos(ip — 0) + rip s'm(ip — 0), (22) 

which implies that 

sin j3 — f sm(f — 0) + rtf> coa(ip ~ 0) . (23) 

Now, using Eqs. 0; an d HSJ and the two- 

dimensional Euclidean line clement in polar coordinates 



The paths followed by light traversing the director config- 
urations corresponding to the defects studied are given by 
the solutions of Eq. iJBJ These equations were obtained in 



dl 2 = dr 2 + r 2 



we get 



ds' 



= { 

K 

{2(, 



2 2 

n„ cos 



1)0- 



sm 

2 



[(*" 
2 [{k-l) 



n 2 ) 2 sin[(fe 



2 • 2 

n„ sm 



? 2 

n~ cos 



[(*- 
c] cos[(A; — 



(24) 



l)<j> + c]}dr 2 
<j> + c]}r 2 d(f> 2 
1)0 + c]} rdrdcj). 

(25) 



For the metric l|25|l the only nonzero Riemann curva- 
ture tensor components are given by 



R rtj}r4> = -k{k - l)(n 2 e - nl) cos{2[(fc - 1) 



and permutations R<f, r ft 
The curvature scalar is 



-R, 



<firr<fi 



c]}. (26) 



R = 



2k{k-l){nl-nl) 



cos{2[(fc- 1)0 + c]}. (27) 



By inspection of Eqs. (|25|l . I|2t)|l and (|27|l we see that, 
in the isotropic limit, n — n e , as in the hedgehog case, we 
recover Euclidean space. Also, for n ^ n e k — 1 and c = 
and c = |we reproduce the result of > which associates 
the effective geometry of the k = 1 defects with that of a 
cosmic string. A similar result was found for disclinations 
in elastic solids in the geometric theory of defects approach 
|18| . In particular, the k — 1 and c = case coincides also 
with the 8 — tt/2 section of the hedgehog effective metric 
(Eq. jnj 



ds = n n dr 



n 2 r 2 c 



(28) 



This implies that light propagating in the presence of the 
hedgehog, on the surface 9 — tt/2, has the same behavior 
as light propagating on the surface z = const of the k = 1, 
c = disclination 





Fig. 5. Director field and light paths for the 9 = tt/2 section 
of the hedgehog (same as k — 1, c = disclination) and k = 1, 
c = J disclination, respectively. 





Fig. 6. Director field and light paths for the k — — 1, 
disclination. 





Fig. 7. Director field and light paths (vertical) for 
c = ? and k = — i c = disclinations, respectively. 
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Fig. 8. Director field and light paths (skew) k — 1 
k — — i, c = disclinations, respectively. 
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Fig. 9. Director field and light paths (horizontal) k — = , < 
and fc = — 1, c = disclinations, respectively. 
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their final form after the computation of the Christoffel 
symbols (Eq. Q) obtained from the metrics (Eqs. Ijl7|l 
and l |25jl and solved numerically using the Runge-Kutta 
method. For the numerical calculation we used the same 
values as reference jSJ: n = 1.65 and n e = 1.94 corre- 
sponding to the nematic Merck Phase V and look only at 
the extraordinary wave. In fact, the exact values are not 
important because the behavior obtained will be the same 
as far as n < n e . 

We consider the propagation of a iniatially parallel 
light beam through the defect region in order to get a 
feeling for the bending of the rays by the defects. Figs. 
5-9 show assorted lensing behavior from the defects stud- 
ied. In particular, Fig. 5 shows a clear lensing behavior 
for both k — 1 defects, the c = defect acting as a di- 
verging lens and the c = ir/2 defect acting as a converging 
lens. The converging lens behavior is similar to the one 
presented by the cosmic string and by the irrotational 
hydrodynamical vortex [HJ. The diverging lens behavior 
is similar to the one presented by a cosmic string with 
negative mass density [j]]. Like the vortex and differently 
from the cosmic string, the lensing effect here appears as 
consequence of the induced geometry, the effective geom- 
etry perceived by the light rays crossing the medium. The 
k = — 1 and k = ±1/2 defects, being asymmetric, have 
different behavior depending on the incidence angle of the 
light, as shown in Figs. 6-9. 

The overall behavior observed in Figs. 5 - 9 is consis- 
tent with the the variation of the refractive index with 
the angle j3 between the ray and the optical axis of the 
molecule. Notice that the optical length between two given 
points, which should be minimal, will be smaller for smaller 
N r (see Eqs. © and (JTJ). Therefore, the ray "chooses" to 
travel along the paths with smaller N r , which in this case 
is parallel or antiparallel to the molecule (see Fig. 10). 



mm 



4 Concluding Remarks 

We study the propagation of light in a uniaxial liquid crys- 
tal for the specific molecular configurations corresponding 
to selected topological defects. We use a differential geo- 
metric approach, where the light rays are considered as 
geodesies in a generic orientation-dependent space. Natu- 
rally, the liquid crystal is in flat Euclidean space, only the 
bending of the extraordinary rays is interpreted as due to 
geometry. In other words, we ask ourselves what is the ge- 
ometry that gives the same paths as the ones followed by 
extraordinary light rays in nematics with topological de- 
fects. As a result we find a striking similarity with cosmic 
defects: the effective geometry associated to the hedge- 
hog defect is that of the global monopole and that of the 
k = 1 disclinations is that of the cosmic string. The light 
paths near the defects are geodesies obtained from the 
geometry and were computed numerically showing an ob- 
vious lensing effect, in agreement with the earlier result of 
Grandjean [7|. This similarity suggests that experiments 
to study optical properties of cosmic defects can be done 
with liquid crystals in the laboratory. 

Analogue gravitational systems in condensed matter 
allow for laboratory tests of some hypotheses otherwise of 
very difficult realization with present technology. Among 
these is gravitational lensing, an important research area 
in astronomy and astrophysics, which may benefit from 
the analogy shown in this article between cosmic defects 
and defects in liquid crystals. A major difficulty in the 
research on gravitational lensing is the identification of 
the nature of a lens candidate. The possibility of simu- 
lating gravitational lensing behavior in the laboratory ap- 
pears as a tool that might help in the characterization 
of known gravitational lenses like, for example, the CSL- 
1 (Capodimonte-Sternberg-Lens Candidate, No.l), an ex- 
tragalactic double source believed to be a cosmic string 
19 . Although asymmetric cosmic strings are not yet known, 
a comparative study of the known multiple images formed 
by the many known gravitational lenses and the images 
formed by the k = ±| defects may indicate the presence 
of such objects. 



5 Appendix: Finslerian Approach 




Fig. 10. Refractive index N r as a function of the angle /3 
between the ray and the molecule's optical axis, from Eq. Q. 
(3 is given in radians. 



In Finsler geometry , the line element depends not only 
on the position but also on the infinitesimal elements of 
direction dx % , which can be written as dx l = x l dt, where dt 
is a parameter along some curve. If we think of t as time, x l 
is the velocity of a particle moving along the curve which 
is therefore tangent to the curve. The arc length element 
is given by 

ds = F(x,x)dt. (29) 

where F(x,x) is the Finslerian function. The geodesic 
equation in Finslerian geometry is 



r 



dx j dx k d(logF) dx % 



dt dt 



dt dt 



(30) 
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In the Joets and Ribotta approach to geometric optics 
in anisotropic media [H] the 

F(x, x) — N r (31) 

identification is done. Equation H3()(l can be simplified by 
noticing that F will be constant (i.e. independent of t) 
and equal to 1 if t is taken to be the Finslerian arc length 
s along the geodesic, which follows immediately by mak- 
ing dt — ds in Eq. I|29|) . We then recover equation JHJl of 
Riemannian geometry. The Christoffel symbols are com- 
puted in the same way as in Eq. (JJJ and the Finslerian 
metric tensor components, 

1 d*F* 

913 ~ 2 dx^dxJ ' [ ' 

turn out to be exactly the same ones computed with the 
help of Riemannian geometry for the cases studied. This 
justifies the use of Riemannian geometry for the calcula- 
tions presented in the main body of this article. 
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